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Abstract. This article shows that the Neumann dynamical sys- 
tem is described well in terms of the Weierestrass hyperelliptic al 
functions. 

1. Introduction 

The Neumann dynamical system is a well-known integrable nonlinear 
dynamical system, whose Lagrangian for (q, q) e IR 29+2 is given by, 

9+1 9+1 

(1-1) L = \^.--Y,^l 

i=l i=l 

with a holonomic constraint, 

9+1 

(1.2) <%)=0, $(g):=^g 2 -l. 

i=l 

This is studied well in frameworks of the dynamical system |Moj . of 
the symplectic geometry |GSj . of the algebraic geometry |Muj . of the 
representation of the infinite Lie algebra [AFlPl IS]. 



Mumford gives explicit expressions of the Neumann system in terms 
of hyperelliptic functions based upon classical and modern hyperelliptic 
function theories. This article gives more explicit expressions of the 
Neumann system using Weierstrass hyperelliptic al functions. 

In the case of elliptic functions theory, Weierstrass p functions and 
Jacobi sn, cn, dn functions play important roles in the theory even 
though they are expressed by the 9 functions and all relations among 
them are rewritten by the 6 functions. The expressions of Weierstrass 
p functions and Jacobi sn, cn, dn functions make the theory of elliptic 
functions fruitful and reveal the essentials of elliptic functions. 

Unfortunately in the case of higher genus case, such studies are not 
enough though Klein and Weierstrass defined hyperelliptic versions of 
these p functions [KT| and sn, cn, dn functions jW]. Thus several 
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authors devote themselves to reinterpretations of the modern theory 
of hyperelliptic functions in terms of these functions and developing 
studies of these functions as special functions [|BEL[ IMaj and their 



references]. In this article, we also proceed with such a project. We 
will show that the Weierstrass al functions give natural descriptions of 
the Neumann dynamical system: As in theorem 4.1, the configuration 
q % of z-th particle (or coordinate) is directly given by the al function, 

q\t) = ak(t), 

Here aU(t) are defined in Definition 3.1, which was originally defined by 
Weierstrass as a generalization of Jacobi sn, cn, and dn functions over 
a elliptic curve to that over a hyperelliptic curve. As Jacobi sn, cn, dn 
functions are associated with several nonlinear phenomena and these 
relations enable us to recognize the essentials of the phenomena [Tj , we 
expect that this expression also plays a role in hyperelliptic function 
case. In fact the description in terms of the al functions makes several 
properties of the Neumann system rather simple. For examples, an 
essential property of the Neumann system (? 4 (t)) 2 = 1 is interprited 
as a hyperelliptic version of sn 2 (w) + cn 2 (w) = 1. Its hamiltonian is 
given as a manifestly constant quantity in Theorem 4.1 (3). Due to the 
description, proofs in this article basically need only primitive residual 
computations. 

We will give our plan of this article. §2 gives a short review of the 
Neumann system. In §3, we introduce the hyperelliptic al functions 
and hyperelliptic p functions. There we also give a short review of 
their basic properties following |Ba| IBEl"| IWj. §4 is our main section, 
where we give our main theorem. There q[ function naturally describes 
the Neumann system. 

We thank Professor Emma Previato for bringing my attentions up on 
the Neumann system and Professor Yoshihiro Onish for his continual 
supports on the studies. 



2. Neumann System 

We will shortly review the Neumann system (q, q) G M 29+2 whose 
Lagrangian and constraint condition are given (jl.lj) and (jl.2|) in Intro- 
duction. The constraint (jl.2|) means <&(q) = 0, 

9+1 

(2.1) = 0. 

i=i 
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The canonical momentum pi to ft is given as 

dL . 

Pi = ^- = q- 

oqi 

Purely kinematic investigations lead the following proposition [M 
Proposition 2.1. The hamiltonian of this system is given by 

9+1 9+1 



(2.2) ff^Etf + ^E**?' 



2^" 2 

i=i i=i 

and i/ie hamiltonian vector field is given by 

(2.3) = E^|- - E a ^|- + (E^ 2 - #0 E*|:- 

TTie equation of motion is given by 

(2.4) ft = ft, ft = -(2L + aj)ft. 



3. HYPERELLIPTIC FUNCTIONS 

In this article, we will consider a hyperelliptic curve C g given by an 
affine equation |Mu| IDRVWj , 

y 2 = f(x), f{x)=A{x)Q{x), 

A(x) := (x - a 1 )(x - a 2 ) • ■ ■ (x — a g+1 ), 
Q(x) := (x - ci)(x - c 2 ) ■ ■ ■ (x - c g ), 

where Oj's and Cj's are complex numbers. Let hi :— a,i (i = 1, • • • , g + 1) 
and b g+i+1 := c 4 (i = 1, • • • , g). 

In this article, we deal with {x\, x 2 , • • • , x g ) belonging to g symmetric 
product Sym 9 (C g ) of C g . 

Let us introduce the canonical coordinate u := {u\, • • • , u g ) in the 
Jacobian J g related to C g BEE], 

9 r(x a ,y a ) x i-^d x 



Ui := E 

a=l 



2y 



Here u_ := (u±, • ■ • , « g -i), u = (it-, u g ). 

Due to Abel theorem jHj, the following proposition holds. 

Proposition 3.1. (ui,u 2 ,--- ,u g ) are linearly independent in C 9 . In 
other words, there are paths in Sym 9 (C 9 ) so that {u g } is equal to C 
with fixing w_ . 
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As Mumford studied the Neumann system using t/VVT-expression 
of the hyperelliptic functions |Muj . we will give U, V and W functions 



(3.1) U(x) =: (x - xi) ■ ■ ■ (x - x g ), 

In this article, we will express the system in terms of the hyperelliptic 
p functions and al functions. Let us introduce these functions, 

Definition 3.1. The hyperelliptic p g i (i = 1, 2, • • • , g) functions ofu's 
are defined by 

g 

(3.4) U{x) = x° + 

i=l 

e.g., p gg := X\ H V x g . 

The Weierstrass ah and aU (i — 1, 2, • • • , g) functions are defined by 

M 



(3.5) al r (u) := 7 r al r (-u), al r (w) := y / U(a r )(u), 



where we set 7 r = 1/ A'(a r ) in this article. We write 
al^u) := ^aL(u), af{u) := ^^(u) 



The sn function is defined by l/sn(u) := \/x — a$j a/ A' (as) and 
sn(u) = l/sn(w + 0), al functions should be recognized as an extension 
of sn function. As sn function has the relations 

k 2 sn 2 (u) + dn 2 (w) = 1, sn 2 (w) + cn 2 (w) = 1. 

the at functions are also has similar relations, which were studied in 
luj as a generalization of Frobenius identity. 



Proposition 3.2 



a 

i=l i=l 



NEUMANN SYSTEM AND HYPERELLIPTIC AL FUNCTIONS 

Proof. The left hand side is given by 

9+1 U(a t ) 1 ^ U(x) 



EU(ai) 1 \ - 



since around the finite ramified point (aj, 0) of the curve C g , we have 
a local parameter t 2 = {x — a») and 

C/(a;) , 2t/(t 2 + Oi)tdt 
res (aii0 )— --dx = res (aii0 )- 



'A(x) (a " u > (t 2 + ai- ai )---t 2 --- {f + ai -a 



9+1, 



Let us consider an integral over a boundary of polygon expression Co 

U(x) 



oiC g , 



dC 



A(x) 



which gives the relation, 



U(x) U(x) 

At oo, a local parameter t of C g is given by x = 1/t 2 : 

U(x) _ ^{l-x l t 2 ) ■■■{!- x g t 2 ) -2^_ 

res OOA(a;)d x-res 00 _^ (i _ ait2) __ _ (i _ a ^ t2) ^ - 2. 

Hence it is proved. Similarly we obtain the relations for alj?' though 
we should evaluate W(x)/x A (x). □ 

There is a natural relation between al function and p gg function 

Proposition 3.3. 

d 2 d ( 29+1 \ 



\j = l,bi^CLi 



Proof. This is directly obtained if we assume the following Lemma 13.11 
andEI3 □ 



We have primitive relations between differentials of al functions and 
UVW expressions: 



Lemma 3.1. 



d 9] M = V{a ^ u \ a&Hu)- V{ai){u) 



ali(w) ' 1 aU(u)A'(ai) 
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\X la) 



Proof. Noting , 

a=i ^ A Xq ) 
1 .9 

we find that - — — U(x) = V(x), which directly gives the relations. □ 



2du c 



Lemma 3.2 

1 d 



2du, 



7(ai) = £/ (ai) I 6< - aj - 2 x a 

V i=l a=l 



C/(ai 



Proof. Here we will check the left hand side. 
—V( ) = V — - ^( a ») 

E 2?/ a g 2y a U(ai) y a 9 2y b U(a.i 



du 



^ U'(x a ) dx a U'(x b )(ai - x b ) 



U'(x a ) dx a U'(x b )(ai - x b ) 



U{ai) it I 



0=1 



1 d ( f{x)U{oi) 



U'(x) dx \U'(x)(ai — x b ) 

f(Xa) 



9 



U'{x a ) 2 {a i -x a y 



2y a 



2/b 



23+1 9 



U'{x a ) U'(x b )(ai - x a )(a,i - x b ) 

9 



2y a 



1=1 



a=l 



U'(x a )(a,i - x a ) 



Here we used the following relations. 



1 9 / f(x)U( ai ) 



U'(x) dx \U'(x)(ai — Xb) 



Tes (xa,y a ) 



and 



23+1 



Y b t - di - 2 Y x a = ^ res (w<0 



U(x) 2 (cii — x) 
fix) 



dx, 



-dx. 



) U{x) 2 {a i -x ) 

i=l a=l a=l \ / \ / 

The third relation is obtained by an evaluation of the integral 

fix 



dc U(x) 2 (di-x 



-dx. 



□ 
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Remark 3.1. The Klein hyperelliptic p function obeys the KdV equa- 

d 

tions IB EL I [Mai. On the other hand, — — logal r is a solution of the 

8Ug 

MKdV equation |Maj . The relation in Proposition 13.31 means so-called 
Miura transformation, 

d \ 2 d 2 

log alj ) + ^ log ah = {C- di), 



du g J du g 



where £ := - (2p gg - k). 



4. Neumann system and hyperelliptic al functions 

This section gives our main theorem as follows. 

Theorem 4.1. Suppose that configurations of(xi, ■ ■ ■ ,x g ) G Sym 9 (C ff ) 
are given so that (alj) belongs to M. 9+1 , u 9 6l fixing w_ G 

(1) ali obey the Neumann system, i.e., 

(4.1) 9i (f) = (*(«_,*), ft = ot ] («_,t), 

where the time t of the system is identified with u g and thus the 
Hamiltonian vector field is given by 



dt du g 



(2) The hamiltonian \2.2]) and the lagrangian M.l\) are given by 

1 ( g+l 9 \ if 29+1 \ 

H = 2 ^ Gi ~ Ca ' L= 2\ 2p " ~ ^ bi I ' 

\i=l a=l / \ 1=1/ 

(3) The conserved quantities are q (i = 1, • • • , g) and 

9+1 3+1 / . • \2 

mi := % + > > , (« = 1, • • • , g + 1), 

< • < • n ■ — n ■ 



i=l j=l,^i 

which obey relations, 

3+1 9+1 

iii, = — — -. y]mj = 1, ^ ",///, //. 



These relations were essentially proved in |Muj using UVW expres- 
sion without al functions. However from the viewpoint of studies of 
special functions, we will show them directly using nature of al func- 
tions. 
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Proof. Assumptions are asserted bv Proposition lTTl (1): Due to Propo- 
sition I3~2l ai r 's obviously obeys the constraint condition $(oQ = (jl.2j) 
and $(at) = (j2.1j) by differentiating the both sides of the identity in 
u g . We should check whether they obey the equation of motion (|2.4|) . 
which are proved in Proposition 13.31 if we assume the form of the La- 
grangian L in (2). (2) is directly obtained by using the relations in 
Lemma [4. 11 Finally (3) is proved in Remark 14.31 □ 

Remark 4.1. (1) The equation of motion (|2.4j) is directly related 
to Proposition ^. 31 which is connected with the Miura transfor- 
mation. Further the constraint (|1.2|) satisfies due to the identity 
of at function as mentioned in Proposition 13.21 These exhibits 
essentials of al functions. Hence the Neumann system should 
be expressed by the al function as some dynamical systems are 
expressed by Jacobi sn, cn, dn functions jTj. 

(2) We remark that the hamiltonian depends only upon Oj's and 
Cj's which determines the hyperelliptic curve C g . Thus it is 
manifest that it is invariant for the time u g development of the 
system. 

(3) There are 2g degrees of freedom as a kinematic system because 
the constraints $ and $ reduce (2g + 2) ones to 2g ones. The 
independent conserved quantities rrii are g — g + 1 — 1; " — 1" 
comes from Y2 m i = 1- Since the sum of rrii gives hamiltonian 
H, H is not linearly independent conserved quantities. Since 
there are other g conserved quantities Cj but their sum gives the 
hamiltonian ^ m «5 the dimensional of independent q is g — 1. 
However Yl 9 i=iQi/ a i — 1 compensates the lacking one. Hence 
the degrees of freedom of this system is equal to number of the 
conserved quantities. 

(4) By the definition of q's, q depends upon the initial condition of 
the Neumann system whereas is fixed as coupling constants 
of the Neumann system. Thus S g := {C g : y 2 = A(x)Q(x) \ 
Ci,c 2 , • • -c g G C} corresponds to the solution space M g of the 
Neumann system if u g E R and (at, ai [g] ) G R 29+2 . The S g is a 
subspace of the moduli M. g of hyperelliptic curves of genus g. 

Let us give a lemma and remarks as follows, which are parts of the 
proofs of the theorem. 



9+1 9 

Lemma 4.1. (1) ^[at| s] ( M )] 2 = p 9g (u) - ^c a . 

i=l a=l 
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3+1 

(2) S ^2a i a{ i (u) 2 = y^Qj - p gg (u). 
1=1 

/ V(x) 2 

Proof. 1) Due to Lemma ETT1 we deal with <b dx = giving 

Jdc a U{x)A{x) 

^ ^(a^) 2 V(x) 2 V(x) 2 

2 U( ai )A'( ai ) + a J^ ± reS ^'^) rfx+res °° ^(x) ^ = °' 

Whereas the third term vanishes, each element in the second term is 
given by 

reS(Xa ' ±ya) ?7(x)A(a;) tf'(x a )' 
/ Q(x) 

Further we also evaluate an integral, d> — ——dx = 0. The integrand 

Jdc U{x) 

has singularities at (x a , ±y a ) and the infinity. Similar consideration 
leads us to the identities 

£ §H = ^ 1 + --- + c g )-2(x 1 + ... + x g ). 

a=l,e=± U ^ 

Due to these relations, we have the finial equal. 

/ U(x) 

2) Next we will consider an integral, <p x dx = 0. A residual 

EU (di) U(x) . . 

ctj— — — - = —res 00 x———dx. lhe infinity term 
, =i A'{ai) A{x) 

gives 2((xi + ■ • • x g ) — (ai H %+i))- Hence we also have the relation 

in (2). ' ' □ 

dXd ^Va 

Remark 4.2. Using the fact — — = — y— —, we obtain anther form of 

du g U'[x a ) 

Lemma\£l$Mm, Y}^f = E ^^tp where g(x) a , b : = 

i=l a, 6=1 " ^ 

— — ; ■ — — — — — - whose off-diagonal part does not vanish 

[ai- x a ){ai- x b )A'{ai) 

for the case genus g > 2 in general. 

Remark 4.3. (Proof of Theorem 4-1 3) Here we will give the conserved 
quantities of the Neumann system. Let us consider, 

9+1 g+l , . . s 2 

mi{x) = Qi + \ \ J - J - . 

^— ' / — ' X — CLj 

1=1 j=l,^i J 
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Then we have identities 



f(x) = U(x)W(x) - V(xf _ mi (x) 
~A{xf = A{x) 2 ~ ^^a? 

m A x ) 2 . ~~ W 



mi = res ai = g. 



x - a, ; f—f a« - a 

1=1 J = l,7±l 



J 



The direct computation gives the relations in Theorem when we deal 

with the integrals of differentials ^ \ dx, ®. . dx. 

A(x) A(x) 
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